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Abstract 

A  Reynolds-average  Navier  Stokes  method  developed  for  rotorcraft  type  of  flow  problems  is  applied  for 
predicting  the  forces  and  moments  of  forebody  models  at  high-angle-of-attack  dynamic  conditions  and  for 
providing  insight  into  the  flow  characteristics  at  these  conditions.  Wind-tunnel  results  from  rotary  testing 
on  generic  forebody  models  conducted  by  NASA  Langley  and  DERA  are  used  for  comparison.  This 
paper  focuses  on  the  steady-state  flow  problem. 

Introduction 

The  stability  and  control  at  post-stall  maneuvering  flight  conditions  are  of  major  concern  in  the  design  of 
advanced  agile  airplanes.  At  these  flight  conditions  the  flows  about  this  type  of  aircraft  tend  to  be 
unsteady  and  the  aerodynamic  forces  and  moments  show  strong  nonlinear  variations  with  rotation  rates. 
During  the  development  process  of  these  airplanes,  numerous  hours  are  spent  testing  scaled  models  in 
ground  facilities  with  complex  special  purpose  equipment  to  determine  the  dynamic  stability  and  control 
characteristics.  Much  of  this  testing  is  conducted  using  sub-scale  models  of  simplified  geometry  at  low 
Reynolds  numbers.  In  addition  to  this  largely  experimental  approach,  potential-flow  type  of  prediction 
methods  including  empirical  viscous  corrections  are  being  used  to  predict  loads  on  forebodies  under 
separated-flow  conditions.  However,  the  sensitivity  of  the  forces  and  moments  to  geometric  detail  and 
Reynolds  number  makes  these  (semi-)empirical  approaches  to  predict  the  dynamic  characteristics  at  full- 
scale  high-angle-of-attack  conditions  a  challenging  task. 

In  the  mid-1990s,  rotary  balance  experiments  were  conducted  on  square  and  circular  ogive  forebodies  at 
angles  of  attack  of  60  and  90  degrees  over  a  range  of  Reynolds  numbers  (based  on  the  body 
diameter/width,  D,  and  freestream  velocity,  Uoo)  from  0.08  x  106  to  2.25  x  106  (Refs.  1  and  2).  Ligure  1 
depicts  the  two  forebody  models  with  square  and  circular  cross  sectional  shape.  The  purpose  of  these 
experiments  was  to  determine  the  effects  of  Reynolds  number,  angular  velocity,  and  geometric  shape  on 
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the  aerodynamic  characteristics  of  the  forebodies.  These  tests  were  unique  in  so  far  that  this  was  the  first 
time  that  surface  pressure  distributions  were  measured  under  dynamic  conditions  in  a  pressurized  wind 
tunnel.  A  second  objective  of  these  experiments  was  to  provide  a  database  for  the  development/validation 
of  high-angle-of-attack  computational  methods.  This  database  forms  the  basis  for  the  computational 
investigation. 

The  numerical  simulation  of  rotary  aerodynamics  is  only  in  its  infancy.  Inviscid  methods  and 
viscous/inviscid  interaction  methods  tend  to  be  inapplicable  because  of  the  strong  viscous  effects  and  the 
large  regions  of  flow  separation  that  occur  at  high-angle-of-attack  conditions.  Also,  methods  based  on  the 
Euler  equations  only  predict  vortex  formation  originating  from  flow  separation  at  sharp  edges  and  ignore 
the  vortices  that  originate  from  flow  separation  at  smooth  surfaces.  At  present,  the  only  viable 
computational  approach  for  the  prediction  of  aircraft  dynamics  at  high  angles  of  attack  is  based  on  the 
time-dependent  three-dimensional  Reynolds-averaged  Navier-Stokes  (RaNS)  equations.  Unfortunately 
much  of  the  developments  in  computational  fluid  dynamics  (CFD)  have  concentrated  on  static  problems, 
such  as  the  computation  of  lift  and  drag  of  aircraft  at  stationary  conditions.  Relatively  few  researchers 
have  explored  the  application  of  three-dimensional  CFD  methods  to  aircraft  dynamic  problems. 

A  few  studies  provide  approaches  to  deal  with  body  motion  in  CFD  simulations  including  Kandil  & 
Chuang  (Ref.  3),  Weinacht,  Sturek  et  al  (Refs.  4-7),  Chaderjian  &  Schiff  (Ref.  8),  Ahmad  &  Duque  (Ref. 
9),  Limache  &  Cliff  (Ref.  10),  and  Park  &  Green  (Ref.  11).  Kandil  &  Chuang  (Ref.  3)  studied  the  flow 
about  an  oscillating  delta  wing  at  supersonic  conditions  using  the  RaNS  equations  in  conical  form. 
Weinacht,  Sturek  et  al  (Refs.  4-7)  analyzed  the  damping  characteristics  of  various  motions  (including 
coning  motion)  of  projectiles  at  supersonic  conditions  using  a  parabolized  Navier-Stokes  method. 
Chaderjian  &  Schiff  (Ref.  8)  analyzed  the  wing-rock  problem  of  delta  wings  at  high  angle  of  attack  using 
a  single-grid  technique  whereas  Ahmad  &  Duque  (Ref.  9)  analyzed  the  flow  around  a  rotating  helicopter 
rotor  using  a  more  general  overset  grid  technique.  Through  various  enhancements  of  the  RaNS  solver 
OVERFLOW,  the  latter  methodology  has  been  developed  further  by  the  Army  /NASA  rotorcraft  group  at 
Ames  Research  Center.  These  enhancements  have  made  it  feasible  to  efficiently  simulate  the  rotary  flows 
over  forebodies  without  the  need  for  extensive  modifications  of  flow  solvers  developed  for  static 
problems.  Limache  &  Cliff  (Ref.  10)  analyzed  rotary  stability  derivatives  using  an  ingenious  sensitivity 
equation  method  and  applied  this  scheme  to  solve  the  pitch-rate  derivatives  for  airfoils  at  inviscid  flow 
conditions.  Park  &  Green  (Ref.  11)  modified  the  general  purpose  RaNS  solver  CFL3D  to  compute 
stability  derivatives  for  a  various  configurations  including  a  complete  airplane. 

In  the  next  section  the  forebody  models  as  well  as  the  experimental  setup  used  to  generate  the  database 
are  discussed,  followed  by  a  brief  discussion  of  the  rotorcraft  version  of  OVERFLOW  and  the 
computational  grids  on  which  the  solutions  were  obtained.  Computational  results  are  presented  for  the 
circular  and  square  ogives  at  steady  rotary  conditions  and  these  results  are  compared  against  the 
experimental  findings. 


Geometry  of  Forebod v  Models  and  Test  Setup 

The  rotary  flow  experiments  were  part  of  a  cooperative  effort  between  NASA  Langley  and  DERA 
(Defense  Evaluation  &  Research  Agency)  and  were  conducted  in  the  United  Kingdom  in  the  atmospheric 
13  ft  x  9  ft  wind  tunnel  at  Bedford  and  in  the  pressurized  8  ft  x  6  ft  wind  tunnel  at  Famborough,  utilizing 
the  DERA  rotary  balance  apparatus.  These  wind  tunnels  provided  a  freestream  Mach  number  range  from 
0.024  to  0.21  and  a  Reynolds  number  range  from  0.15  to  4.5  million  per  foot  for  these  tests.  The  sketch 
in  Fig.  2  depicts  the  test  setup  as  it  was  numerically  simulated  in  the  present  investigation. 

Detailed  descriptions  of  the  forebody  models  are  presented  by  Pauley  et  al  (Ref.  1)  and  Dunham  (Ref.  2) 
and  sketches  of  the  models  are  shown  in  Fig.  1.  As  shown  in  Fig.  1,  one  set  of  model  configurations  has  a 
square  cross  sectional  shape  with  rounded  corners  whereas  the  second  set  of  models  has  a  circular  cross 
sectional  shape.  The  detachable  aftbodies  have  a  hemispherical  shape.  The  detachable  forebodies  have 
either  a  hemispherical  or  a  2.0D  ogive  shape.  Although  both  the  square  and  the  circular  configurations 
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were  wind-tunnel  tested  with  the  hemispherical  as  well  as  the  ogive  forebodies,  only  the  configurations 
with  the  ogive  forebody  shape  are  considered  in  this  computational  study.  Detailed  surface  measurements 
were  conducted  for  both  ogive  forebody  models  and  were  used  to  create  the  surface  grids.  Both 
forebodies  turn  out  to  be  slightly  asymmetric. 

The  ogive  models  have  six  rows  of  pressure  taps  in  the  forebody  and  two  rows  of  pressure  taps  in  the  aft 
section  for  a  total  of  254  taps.  Each  model  has  a  length,  b,  of  36  in.  and  a  maximum  diameter/width,  D, 
of  6  in.  The  forebody  has  a  length  of  12  in.  and  the  hemispherical  aftbody  a  length  of  3  in.  Each  model  is 
mounted  to  the  internal  force-and-moment  strain  gage  balance  which  in  turn  is  connected  to  the  rotary 
apparatus  via  a  sting  located  18  in.  from  the  nose.  Note  that  neither  the  sting  nor  any  other  rotary-balance 
equipment  is  included  in  the  computational  model  at  this  point  in  the  study. 

Flow  Solver 

The  Reynolds-average  Navier  Stokes  (RaNS)  method  used  to  solve  the  flows  about  forebodies  at  high- 
angle-of-attack  dynamic  conditions  is  based  on  version  1.8  of  OVERFLOW  (Ref.  12).  OVERFLOW  has 
been  developed  by  Buning  et  al  for  problems  involving  fixed-wing  aircraft  (Ref.  12).  It  solves  the 
governing  flow  equations  on  structured  overset  grids.  The  reasons  for  selecting  this  code  are  multifold. 
First,  the  code  is  well  documented  and  supported  by  Buning  and  his  colleagues.  Second,  the  code  uses 
the  overset  or  Chimera  grid  approach  to  facilitate  the  solution  of  flows  involving  complex  geometries. 
Third,  researchers  of  the  Army/NASA  Rotorcraft  Division  at  the  Ames  Research  Center  have  enhanced 
the  basic  OVERFLOW  code  by  adding  rotorcraft- specific  capabilities.  The  enhancements  that  facilitate 
the  solution  of  rotorcraft  type  of  flow  problems  include  rotational  rigid-blade  motion,  dynamic 
interpolation  of  solution  quantities  between  rotating  and  non-rotating  grids  (Ref.  13),  and  a  source-term 
formulation  to  allow  steady  flow  solutions  of  the  rotor-in-hover  problem  (Refs.  14-16).  In  the  present 
study,  particularly  this  last  feature  turned  out  to  be  very  useful  and  resulted  in  significant  savings  in 
computing  time. 


The  resulting  Navier-Stokes  equations  for  a  rotating  body  in  a  compressible  flow  can  be  written  as 
follows  (Ref.  16): 


^+E(Q)=V(Q)+R(Q) 


(1) 


where  Q  =  [p,pu,pv,pw,e]T  is  the  solution  vector  of  conserved  variables  for  the  mass,  momentum,  and 

3E 

energy  equations.  Also  in  Eq.  (1),  the  variable  e(q)=-— L  represents  the  convective  operator, 


dx- 


dV; 


V(Q)=  1 1  represents  the  thin-layer  viscous  operator,  and  R(Q)  =  [0,pvfl,-pufl,0,0]Trepresents  the  source 


term  to  allow  solution  of  the  rotary  flow  problem  using  the  steady  formulation  of  the  equations.  Note  that 
this  formulation  fixes  rotation  at  constant  angular  velocity,  f 2 ,  to  that  about  the  z-axis  (Fig.  2). 


OVERFLOW  provides  a  wide  choice  of  algorithms  to  solve  the  governing  equations.  The  results 
presented  in  this  paper  are  conducted  fully  turbulent  using  the  Baldwin-Lomax  turbulence  model  (Ref. 
17)  or  the  Spalart-Allmaras  turbulence  model  (Ref.  18).  In  all  cases,  central  differencing  is  used  on  the 
right-hand- side  of  the  equations.  The  LU-SGS  scheme  is  used  on  the  left-hand  side.  This  scheme  was 
selected  because  it  had  been  previously  tested  with  the  source  term  in  place  for  steady  flow  conditions 
(Ref.  16).  The  second-  and  fourth-order  dissipation  terms  are  set  at  the  default  settings  for  subsonic 
viscous  flow  problems;  0.0  and  0.04,  respectively. 

Computational  Grids 

A  total  of  22  axial  stations  supplied  the  information  for  the  surface  grid  generation  for  each  configuration. 
The  surface  grids  for  the  two  configurations  were  generated  with  the  software  package  OVERGRID 
(Refs.  18,  19).  OVERGRID  is  a  graphical  user  interface  for  communicating  with  a  wide  variety  of  grid 
generation  modules.  Although  it  is  tailored  for  overset  grids,  only  single  grids  were  used  to  simulate  the 


(SYA)  39-4 


flows  over  the  circular  and  square  forebody  models.  Each  surface  grid  consists  of  181  equidistant  points 
circumferentially  at  a  total  of  130  and  131  axial  stations  for  the  circular  and  square  ogive,  respectively. 

The  volume  grids  for  the  ogive  forebody  models  were  also  generated  using  OVERGRID.  For  both 
models  the  outer  boundary  of  the  grid  was  specified  at  50  body  diameters.  The  first  four  grid  points  off 
the  surface  have  a  constant  spacing  of  0.00033D.  Starting  at  grid  point  number  five,  the  grid  is  stretched 
at  a  constant  factor  of  approximately  1.25  resulting  in  a  total  number  of  grid  points  of  54.  Close-ups  of 
the  grids  for  the  two  models  are  shown  in  Fig.  3.  These  grids  were  initially  developed  for  ReD  =  2.0  x 
105.  Extensive  validation  studies  for  laminar  flow  over  a  circular  cylinder  resulted  in  grid  independent 
solutions  and  excellent  agreement  with  experimental  data  in  terms  of  Strouhal  number  and  average  drag 
coefficient.  Further  grid  refinement  may  be  required  for  the  higher  Reynolds  number  turbulent  flow  cases 
considered  here. 


Results 

Simulations  were  conducted  for  both  the  circular  ogive  and  the  square  ogive.  The  aerodynamic 
coefficients  are  based  on  freestream  dynamic  pressure  ypJVloo2^,  area  Db  =  216  in2  for  the  forces  and 
moments,  and  length  D  for  the  moments. 

The  circular  ogive  at  a  =  60°,  ReD  =  2.08  x  106,  =  0.21,  £2b/(2Uoo)  =  -0.20  was  analyzed  using  the 

steady  formulation  of  the  Reynolds-averaged  Navier  Stokes  equations  and  fully  turbulent  flow.  The 
turbulence  model  used  was  the  algebraic  Baldwin-Lomax  model  (Ref.  17).  Figures  4  and  5  depict  the 
convergence  characteristics  in  terms  of  the  L2-norm,  the  normal  force  coefficient,  the  side  force 
coefficient,  and  the  yawing  moment  coefficient.  After  9000  iteration  steps,  the  residual  in  terms  of  the 
F2-norm  is  reduced  by  approximately  three-orders  of  magnitude  and  the  computed  force  and  moment 
coefficients  begin  to  level  off.  The  normal  force  coefficient  converges  to  a  value  of  -0.33,  the  side  force 
coefficient  to  a  value  of  0.18,  and  the  yawing  moment  coefficient  to  a  value  of  0.027.  For  this  case,  the 
side  force  and  moment  results  based  on  the  measurements  in  the  pressurized  8  ft  x  6  ft  wind  tunnel  at 
Farnborough  are  reported  in  Ref.  1.  According  to  Ref.  1,  the  experimental  results  indicate  a  side  force 
coefficient  of  0.18  and  a  yawing  moment  coefficient  of  0.012;  i.e.,  excellent  agreement  is  obtained  for  the 
side  force  whereas  the  yawing  moment  is  significantly  overpredicted.  In  Figs.  6  and  7  the  measured  and 
predicted  forebody  pressure  distributions  are  compared.  No  pressures  were  measured  near  the  nose  tip 
and  hence  the  pressures  in  this  region  are  left  undefined.  The  experimental  surface  pressures  were 
determined  using  as  a  reference  the  pressure  in  a  container  located  inside  the  ogives.  During  the  wind 
tunnel  tests,  problems  were  encountered  correlating  this  reference  pressure  to  the  freestream  pressure. 
This  may  explain  the  consistently  high  experimental  pressures  and  the  resulting  discrepancy  between 
experimental  and  computed  pressure  distributions.  The  shift  in  the  experimental  results  is  particularly 
apparent  in  the  flow  attachment  region  where  Cp  values  in  excess  of  unity  are  indicated.  Based  on  the  low 
freestream  Mach  number  (M^  =  0.21)  and  the  low  rotational  velocity  (£lb/(2Uoo)  =  -0.20),  pressure 
coefficients  much  greater  than  unity  would  not  be  expected.  However,  overall  the  pressure  distributions 
show  good  agreement  with  both  distributions  depicting  the  skewness  in  the  distribution  under  the 
influence  of  the  angular  velocity  component.  Figure  8  shows  the  axial  vorticity  component  at  four  body 
stations.  Figure  8b  and  8e  depict  the  axial  vorticity  contours  for  the  forebody  and  the  aftbody, 
respectively,  with  both  plots  clearly  showing  the  rotational  effect.  Figures  8c  and  8d  depict  the  vorticity 
distribution  near  the  axis  of  rotation  and  hence,  these  distributions  don't  show  as  much  skewness. 

The  square  ogive  at  nearly  identical  conditions  (a  =  60°,  ReD  =  2.09  x  106,  MK  =  0.21,  £lb/(2Uoo)  =  -0.20) 
was  analyzed  using  the  identical  formulation  and  turbulence  model.  For  this  case,  the  convergence 
characteristics  are  very  similar  to  those  for  the  circular  ogive.  The  normal  force  coefficient  converges  to  a 
value  of  -0.39,  the  side  force  coefficient  to  0.17,  and  the  yawing  moment  coefficient  to  0.29.  Again,  the 
side  force  and  yawing  moment  results  based  on  the  measurements  in  the  pressurized  8  ft  x  6  ft  wind 
tunnel  at  Farnborough  are  reported  in  Ref.  1.  According  to  Ref.  1,  the  experimental  results  indicate  a  side 
force  coefficient  of  0.16  and  a  yawing  moment  coefficient  of  0.036.  Hence  good  agreement  is  obtained 
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for  the  side  force  whereas  the  yawing  moment  is  overpredicted.  In  Figs.  9  and  10  the  measured  and 
predicted  forebody  pressure  distributions  are  compared.  Note  that  the  color  contour  scale  is  as  defined  in 
Fig.  6.  Except  for  the  beforementioned  shift  in  the  experimental  data,  the  pressure  distributions  show 
good  agreement  with  both  distributions  depicting  the  skewness  in  the  distribution  under  the  influence  of 
the  angular  velocity  component.  The  relatively  good  agreement  in  the  forebody  pressure  distributions 
and,  hence,  forebody  forces  &  moments  raises  the  question  as  to  what  is  causing  the  overprediction  of  the 
yawing  moment. 

The  discrepancy  between  the  measured  and  computed  values  for  the  yawing  moment  coefficient 
motivated  us  to  rerun  this  case  with  the  Spalart-Allmaras  (S-A)  one-equation  turbulence  model  (Ref.  18) 
instead  of  the  Baldwin-Lomax  (B-L)  algebraic  model  (Ref.  17).  Murman  &  Chaderjian  (Ref.  21)  report 
slightly  more  accurate  results  when  the  former  model  was  used  to  study  the  separated  flow  about  a 
tangent-ogive  cylinder  at  high  angles  of  attack.  In  Fig.  1 1  the  forebody  pressure  distributions  and  axial 
vorticity  distributions  are  compared.  The  normal  force  coefficient  is  shown  to  converge  to  a  value  of  - 
0.56,  the  side  force  coefficient  to  a  value  of  0.21,  and  the  yawing  moment  coefficient  to  a  value  of  0.31. 
The  increase  in  the  magnitude  of  the  normal  force  coefficient  (-0.39,  B-L  versus  -0.56,  S-A)  is  in  line 
with  the  findings  of  Murman  &  Chaderjian  (Ref.  21)  for  the  ogive  cylinder  at  high  angles  of  attack. 
However,  the  yawing  moment  coefficient  (0.29,  B-L  versus  0.31,  S-A)  is  much  less  affected  by  the  choice 
of  turbulence  model  and  remains  inconsistent  with  the  wind-tunnel  results. 

Discussion 

The  computational  results  for  the  high-Reynolds-number  rotary  flow  cases  presented  in  the  previous 
section  show  good  agreement  with  the  experimental  results  in  terms  of  the  forebody  pressure 
distributions,  normal  force  coefficients,  and  side  force  coefficients.  The  agreement  for  the  yawing 
moment  coefficient  is  less  consistent.  Keeping  in  mind  that  aerodynamic  moments  are  as  or  more 
difficult  to  predict  than  aerodynamic  drag,  it  is  appropriate  to  refer  to  Ref.  22  for  a  list  of  possible  reasons 
for  this  discrepancy.  As  pointed  out  in  Ref.  22,  the  accuracy  of  CLD-based  force  &  moment  predictions 
depends  on  a  large  number  of  factors  including  geometry  representation,  grid  size/density,  flow  solver, 
convergence  level,  transition  prediction,  and  turbulence  modeling. 

Detailed  surface  measurements  of  the  forebodies  were  taken  to  obtain  an  accurate  representation  of  the 
ogive  geometries.  However,  the  sting  (see  Lig.  2)  which  connects  the  model  to  the  rotary  apparatus  may 
have  a  noticable  effect  on  the  flow  development  over  the  mid-  and  aft-body.  At  this  point  the  sting  is  not 
modeled.  This  device  is  expected  to  affect  the  vortical  flow  development  along  the  aft  part  of  the  body. 
This  could  explain  the  apparent  relatively  good  agreement  between  the  predicted  and  measured  pressures 
on  the  forebody  and  at  the  same  time  the  relatively  poor  agreement  in  the  yawing  moment. 

Grid  size/density  is  nearly  always  an  issue  for  3-D  viscous  flow  simulations.  Although  the  speed  and 
memory  size  of  computers  have  increased,  allowing  the  solution  of  larger  flow  problems,  the  need  to 
reduce  turn-around  as  well  as  CPU  time  and,  hence,  cost,  pushes  users  to  obtain  solutions  on  smaller 
grids.  Also  in  this  study  a  tradeoff  was  made  between  grid  size  and  computational  cost  with  the  current 
grid  size  encompassing  approximately  1.3  million  points  with  130  points  axially,  181  points 
circumferentially,  and  54  points  normally.  A  grid  spacing  of  Acf)  =  1°  is  required  in  the  circumferential 
direction  to  more  accurately  capture  the  flow  separation  phenomena  and  ensuing  vortical  flow  along  the 
lee  side  of  the  body.  Also  the  initial  normal  spacing,  Ay0  =  3.3  x  10‘4  D,  although  shown  to  be 
sufficiently  fine  for  flow  problems  with  ReD  =  O(105),  may  have  to  be  smaller  for  ReD  =  O(106)  problems. 

The  flow  solver  and  convergence  levels  are  less  of  an  issue  in  this  case.  The  flow  solver  is  well  validated 
and  the  simulations  were  continued  until  the  forces  and  moments  had  converged. 

Boundary-layer  transition  is  a  formidable  issue  for  high-a  forebody  flows  as  pointed  out  by  Haines  (Ref. 
23).  Even  at  the  highest  Reynolds  number  considered  here,  ReD  =  2.1  x  106,  small  regions  of  laminar 
flow  exist  (particularly  in  the  nose  region)  which  can  have  a  major  impact  on  the  global  flow 
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development.  Currently,  there  are  no  transition  prediction  schemes  for  3-D  RaNS  solvers  and,  hence,  the 
accurate  specification  of  transition  onset  becomes  a  tedious  and  time-consuming  task.  Transition 
prediction  schemes  have  been  successfully  implemented  for  the  solution  of  2-D  problems  (e.g.,  Ref.  24). 
However,  3-D  transition  prediction  is  at  least  an  order  of  magnitude  more  complicated  than  2-D  transition 
prediction. 

The  modeling  of  turbulence  remains  one  of  the  most  problematic  areas  of  CFD-based  force  &  moment 
prediction.  In  this  study,  the  algebraic  Baldwin-Lomax  model  and  the  one-equation  Spalart-Allmaras 
model  have  been  used.  Based  on  the  work  of  Murman  &  Chaderjian  (Ref.  21)  it  was  expected  that  the 
Spalart-Allmaras  model  would  provide  more  accurate  prediction  of  the  forces  &  moments.  However,  this 
work  shows  the  differences  to  be  small.  If  anything,  the  Baldwin-Lomax  based  predictions  agree  slightly 
better  with  the  experimental  results.  This  is  in  agreement  with  the  most  recent  observations  by  Murman 
(Ref.  25)  who  notes  that  the  application  of  one-  and  two-equation  turbulence  models  to  high-a  forebody 
flows  with  strong  vortices  outside  the  boundary  layer  leads  to  excessive  levels  of  turbulent  eddy-viscosity 
in  these  vortices.  These  high  levels  of  eddy  viscosity  lead  to  a  deterioration  of  the  flow  predictions. 
Murman  (Ref.  25)  suggests  the  use  of  a  vortex  filter  that  is  Galilean-invariant  and  computationally 
efficient  to  improve  the  predictions  based  on  the  one-  and  two-equation  turbulence  models. 

Conclusions 

A  Reynolds-averaged  Navier-Stokes  code,  OVERFLOW,  has  been  used  to  predict  the  rotary  flow 
characteristics  of  two  forebody  models.  This  computational  tool  shows  great  promise  for  the  analysis  of 
major  problems  encountered  in  the  application  of  the  rotary-balance  wind-tunnel  test  technique  for 
aircraft  dynamics  including  model  support  interference,  wind-tunnel  wall  interference,  and  Reynolds 
number  effects. 

Future  efforts  will  include  additional  simulations  to  further  assess  the  effect  of  rotational  speed  and  angle 
of  attack  on  the  forebody  flow  development  and  compare  the  results  with  those  in  the  experimental 
database.  Also,  the  effect  of  the  sting  on  the  flow  development  as  well  as  the  computed  forces  and 
moments  will  be  assessed.  Last  but  not  least,  we  would  like  to  evaluate  ways  to  incorporate  a  3-D 
transition  prediction  methodology  in  the  flow  solver. 
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(b)  Square  ogive 

Fig.  3.  Close-up  of  volume  grid  for  forebody 
model. 


Fig.  2  Sketch  depicting  layout  of  rotary  test  of 
forebody  models.  Body-fixed  axis 
system  (X\  Y\  Z')  is  used  in  experiments 
and  axis  system  (X,  Y,  Z)  is  used  in 
numerical  simulations.  Neither  the  sting 
nor  any  other  rotary-balance  equipment  is 
included  in  the  computational  model. 


Fig.  4  Convergence  characteristics  for  circular 
ogive  at  a  =  60°,  ReD  =  2.08  x  106,  = 

0.21,  and  £2b/(2Uoo)  =  -0.20. 


(a)  Circular  ogive 
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Iteration  step 


Fig.  5  Normal  force,  side  force,  and  yawing 

moment  predictions  for  circular  ogive  at 
a  =  60°,  ReD  =  2.08  x  106,  Mx  =  0.21, 
and  Qb/(2UX)  =  -0.20. 


(a)  Computed 


(b)  Measured 

Fig.  6  Measured  and  computed  forebody 

pressure  distributions  for  circular  ogive  at 
a  =  60°,  ReD  =  2.08  x  106,  Mx  =  0.21, 
and  Qb/(2UX)  =  -0.20. 


Pressure  coefficient,  C 
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Fig.  7  Measured  and  computed  forebody  pressure  distributions  for  circular  ogive  at  a  =  60°,  Reo  =  2.08  x 
106,  Mo,  =  0.21,  and  £2b/(2U*)  =  -0.20.  Stations  as  marked  in  Fig.  1. 
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(c)  Centerbody  station  (66/130) 

Fig.  8  Axial  vorticity  contours  for  circular  ogive  at  a  =  60°,  Reo  =  2.08  X  106,  Mx  =  0.21,  and  £2b/(2Uv,) 
=  -0.20. 
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(b)  Measured 

Fig.  9  Measured  and  computed  forebody  pressure  distributions  for  square  ogive  at  a  =  60° 
ReD  =  2.09  x  106,  M.,,  =  0.21,  and  Qb/(2U„)  =  -0.20. 
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Fig.  10  Measured  and  computed  forebody  pressure  distributions  for  square  ogive  at  a  =  60 
Ren  =  2.09  x  106,  Mx  =  0.21,  and  £2b/(2Ux)  =  -0.20.  Stations  as  marked  in  Fig.  1. 
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Fig.  11  Comparison  between  Baldwin-Lomax  (left)  and  Spalart-Allmaras  (right)  solutions  for  square 
ogive  at  a  =  60°,  ReD  =  2.09  x  106,  =  0.21,  and  £lb/(2Uoo)  =  -0.20. 
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Question  by  Dr.  Radespiel:  Does  the  experimental  data  at  high  Reynolds  number  indicate 
boundary  layer  transition  on  the  forebody? 

Answer:  Unfortunately,  the  transition  location  was  not  determined  during  the  rotary  wind  tunnel 
tests  of  the  given  forebody  models.  There  is  some  limited  data  with  boundary-layer  trips  on  the 
forebody.  However,  it  is  not  clear  that  transition  really  occurred  at  the  trips. 
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